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Graphite under uniaxial compression along c-axis: a new parameter relates 
out-of-plane strain to in-plane phonon frequency 
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Stacking graphene sheets forms graphite. Two in-plane vibrational modes of graphite, Ei„ and 
Ejg , are derived from the graphene E 2 g mode, the shifts of which under compression are considered 
as resnlts of the in-plane bond shortening. Values of the Griineisen parameter have been reported to 
qnantify such relation. However, the reason why the shift rates of these three modes with pressure 
differ is unclear. In this work, we introduce new parameters 7/^23 =-0.0131 and =0.0585 to 
qnantify the contribution of out-of-plane strain to the shift of the in-plane vibrational frequencies, 
suggesting that the compression of the 7r-electrons plays a non-negligible part in both graphite and 
graphene under high pressure. 

PACS numbers: 62.50.-p, 63.20.-e, 63.20.dk, 63.22.Rc 


I. INTRODUCTION 


Graphene has been investigated intensively since its 
discovery in 20 04 ^, for its unique mechanical and elec¬ 
trical propertied^is] motivation to study graphene 

under strain has been illustrated in the literature®^. In 
brief, strain can modify the properties of graphene to 
meet specific scientific or technological needs. Therefore, 
strain determination and monitoring are of critical appli¬ 
cation importance and contribute to the characterization 
and understanding of this remarkable material. Strain is 
related to phonon frequencies, which can be directly ob¬ 
tained by experiments. The analysis, rather naturally, 
is two-dimensional. Thus the frequency shifts of the 
graphene £2^ mode are considered as entirely induced 
by in-plane strain. 

Stacking graphene sheets forms graphite. The in-plane 
vibrational modes Ei„ and of graphite are derived 
from the graphene E2g mode, and from the two di¬ 
mensional analysis, all these three modes were expected 


to have the same shift with in-plane strain. 

The problem is that these shift rates are not the same. 
That is what we investiga te he re by density functional 
theory (DFT) calculation j^^b^ l gf graphite under non¬ 
hydrostatic conditions. 

Huang et al. gave the dynamical equat ions of the 
graphene optical phonon modes E2g asSHSlSl 

^ ^ (1) 
p 

where u={ui,U2) is the relative displacement of the two 
carbon atoms in the unit cell, to is the phonon frequency, 
and K is the force constant tensor, which can be ex¬ 
panded in powers of strain as 

KaP = -K^a/3 Y ^ ^ ^aPlm^lm (2) 

Im 

^a^im independent elements because of 

the hexagonal lattice, so Eqs.[^and[^may also be written 
as 
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where A and B are the two independent elements of 
^apim *^0 is the unperturbed frequency. For small 
shifts, — cjg « 2 wo(a; — wg), Thomsen et al. expressed 
the solution to the secular equation of Eq. with the 
Griineisen parameter and shear deformation potential 
{SDP) aP 

7 ipxx 3 ” ^yy) ^ ^P{pxx (^) 

Wg 2 

where 7° = {A + B)/AijJq (the superscript 0 is used to dis¬ 


tinguish from the refined 7 presented later in this paper) 
and SDP = {B- A)/ 2 iol. 

Eq. makes explicit the two-dimensional nature of 
the analysis. The Grueneisen parameter 7 and the 
SDP are the two key paramet ers a nd a number of 
experimentaJ®^^^ and theoreticaP^ papers reported 
work on graphene under strain to define their accurate 
values. The results are shown as 7° and SDP in TA¬ 
BLE |Tj It is worth noticing that Ghandour et al. pointed 
out that the transverse strain eT =0 rather than eT=- 
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vel, where v is the in-plane Poisson’s ratio and El is 
the longitudinal strain, in the case that uniaxial strain is 
applied by flexure of a beam to which a graphene flake 
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adherecP^. 

For graphite, when two adjacent graphene layers are 
considered, we can simply make two copies of Eq. as 
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where C is added to account for the interlayer coupling. 
The longitudinal modes are not coupled with the trans¬ 
verse modes, giving the zero elements. Eq.|^ still applies 
and the weak interlayer coupling is usually neglected. 
Thomsen et aP obtained the corresponding Griineisen 
parameter as 1.59 (presented in TABLE |l]), from the ex¬ 
perimental data of graphite under hydrostatic pressur^i^. 

We are now able to demonstrate the problem more 
explicitly — if the frequency shifts of the in-plane modes 
are induced by in-plane strain alone, as shown in Eq. 
for graphene, the shift rates of E 2 g with in-plane strain 
(therefore the Griineisen parameter) should be the same 
no matter how the strain is applied and for graphite, Ei^ 
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where D and E are the additional two independent pa¬ 
rameters arising from the new degree of freedom along c- 
axis, D accounting for the compression of the 7 r-electrons 
into the sp^ network and E for the coupling between lay¬ 
ers. These are clearer in the solutions to the secular equa¬ 
tion of Eq. as 

^fl) = + {A + B) X Ein-plane + {D + E) X Ezz 

^^ 2 ) “ + {A + B) X Ein-plane + {D + E) X Ezz 

^(3) ~ T {A B') X Ein-plane T {D E) X Ezz 

— ^oiEin') T (A H) X Ein-plane “b E) X Ezz 

where Exx is equated to Eyy for the following three cases 
in this paper and C accounts for the difference of the fre- 
quency of the Ei„ and E^g modes of unstrained graphite. 


( 2 ) 

and E 2 g modes should have the same shift rates as the 
graphene E 2 g, also no matter how the strain is applied 
(hydrostatic or biaxial). This is against the results shown 
in TABLE [T] It also contradicts the results from previous 
studies on graphite that Ei^ shifts faster than E^g^ under 
hydrostatic pressur^i^HIi] 

The different shift rates between Ei„ and Ei 2 g, or 
E 2 g from various experiments could be a consequence of 
the compression of the 7r-electrons changing the in-plane 
bonds. To describe this effect, we suggest refining Eq. 
making a phenomenological extension as 
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The solutions lead to a new parameter 7 ', relating out-of- 
plane strain to its contribution to the shift of the in-plane 
phonon frequencies, added to Eq. [^as 

— 'yi^Exx “b ^yy) T xSEP{Exx ^yy) '~1 ^zz 
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with (7) 
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where 7 = (A -f- B)jAuj^, SDP = {B - A)l2wl, 7 ^^^ = 
(I)pE)l2Ltj^ and 7 ' ( 2 ) = {D—E)I2uiq, for small shifts. In 

_ -^29 

Section |III[ we model uniaxial strain and uniaxial stress 
along the out-of-plane c-axis, and hydrostatic pressure on 
graphite, to quantify all the parameters and then explain 
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TABLE I. The Griineisen parameter and SDP for obtained from various experiments and calculations of graphene and 
graphite are listed. The values in square brackets are the corrections by Ghandour et aZ.^. 7 °(Eq. is the value calculated 
for each case from Eq. where the values of 7 and 7 / are the same for all the cases and the ratio of out-of-plane strain to 
in-plane is obtained with the approximation that graphene and graphite have the same elastic constants. 


Experiments (graphene) 

7° 

SDP 

7 ° (Eq. 


7 

7 / 

€ 22 / i^^xx ~\~ ^yy ) 

uniaxial strain(beam flexure) 

0.69, %.58] 

0.38,%.435] 

1.90 


1.90 

-0.0131 

0 

uniaxial strain(beam flexure) 

1.99,[^1.34] 

0.99,[^1.31] 

1.90 




0 

uniaxial strain(substrate stretch) 

1.5, 


1.90 




0 

hydrostatic pressure 

1.99,® 


1.48 




32.25 

hydrostatic pressure 

2.3,1^ 


1.48 




32.25 

Calculations (graphene) 








uniaxial strain (in-plane) 

1.87,® 

0.92,® 

1.90 




0 

uniaxial strain (in-plane) 

1.83,® 

1.18,® 

1.90 




0 

biaxial strain (in-plane) 

1 . 8 ,® 


1.90 




0 

hydrostatic pressure 

2 . 0 ,® 


1.48 




32.25 

shear strain 


0 . 66 ,® 






Experiments (graphite) 








hyrdrostatic pressure 

1.59, 


1.48 




32.25 


(2) 

the different shifts of the Ei„, modes of graphite and 
the Ei 2 g of graphene under hydrostatic pressure. 


II. METHODS 

Graphite was studied at 0 K using as im¬ 

plemented in the Vienna Ab initio Simulation Pack¬ 
age (VASPj^. The exchange-correlation effects were 
treated within the generalised gradient approxima¬ 
tion (GGA) as parameterized by Perdew, Burke and 
ErnzerhoP^ and the projector augmented-wave method 
pseudopotential^for carbon were used. To reach highly 
accurate results, we used 900 eV plane-wave cut-off en¬ 
ergy, and the reciprocal unit cell was sampled with 
18x18x9 k-mesh. Van der Waals (vdW) effects were in¬ 
cluded using the Grimme methocpSI as implemented in 
the VASP code. The elastic properties were evaluated 
using the stress-strain methocP^. The vibrational fre¬ 
quencies at the Brillouin zone centre, the P point, were 
calculated using the 2x2x2 supercell employing the finite 
displacement method as implemented in the Phonopy 
codeP. 


III. RESULTS AND DISCUSSION 
A. Geometry 

First of all, we obtain the optimized geometry for un¬ 
strained graphite, as the in-plane bond length of a = 
I. 42 A. and the interlayer distance of c = 3.20A. The 


errors relative to the experimental valueP^ are 0.06% 
and 4.6%. The vdW add-on is included, nevertheless 
the interlayer interaction is not so well-described as the 
in-plane covalent bonding. The EDA calculation (with¬ 
out vdW) usually gives a better agreement to the exper¬ 
imental value of the interlayer distance, however this is 
considered to be a coincidence because EDA is a local ap¬ 
proximation which overestimates bonding. To minimize 
the effects of calculating vdW inaccurately , we study 
the bond anharmonicity under compressive strain, where 
the vdW attractive potential plays only a small role com¬ 
pared to the dominant repulsion. The error in the value 
of interlayer distance would not affect the phonon fre¬ 
quency shift rates with compressive strain as much as it 
would under tensile strain. 


B. Hydrostatic compression 

We then model hydrostatic pressure on graphite by 
setting a smaller unit cell volume than the unstrained, 
optimizing the geometry at that certain volume, and cal¬ 
culating the corresponding sp^ bond length, interlayer 

distance, pressure and phonon frequencies. The frequen- 
( 2 ) 

cies of the Ei„ and E^g modes of unstrained graphite are 
1565.2 and 1559.1 cm“^, respectively. The e rrors relative 
to the experiments are 1.4% and 1.3%P2l, yye assume 
that they are linked to the vdW attractive term and so 
they would not affect the shift rates with compressive 
strain. Phonon frequencies are plotted against pressure 
in Fig. as is the standard for presenting experimental 
data. And the pressure, now as a calculation output, is 
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plotted against the input here — the unit cell volume. 
(L) and (T) refer to two orthogonal in-plane vibrations, 
longitudinal and transverse. The frequency difference be¬ 
tween these two under hydrostatic condition is less than 
0.4 cm“^ for both Ei„ and and the shift rates of (L) 
and (T) with pressure are the same in the case of the 
Eiu and . Therefore, here and in the following cal¬ 
culation, we treat the differnce between longitudinal and 
transverse modes as computational error and will study 
the longitudinal modes alone as a representative. Linear 
least square fits give the shift rates with compressive pres¬ 
sure up to 10 GPa at 5.3 and 4.3 cm“^GPa“^ for Ei„ and 
E 2 g modes, respectively. No experimental data for Ei^ 

exists and the shift rates for E^g^ wer d^^ * ^^ *^ 4.1-4.6. In 
the previous theoretical work, Cousins et al. obtained 
4.74 and 4.67 cm“^GPa“^ for Ei„ and moded^, 
while Abbasi-Perez et aZ.l^got 5.0 and 4.3 cm“^GPa“^. 
To summarize, the calculation results are reliable and 
reasonable, with the shift rates with pressure comparable 
to previous work, with clear sublinearity of the frequency 
shift due to the pressure dependence of the elastic con¬ 
stant Css, and with the two in-plane modes degenerate 
when the graphene layers are pulled apart. However, the 
problems are again the different shift rates for Ei^ and 
^ 2 g with pressure, and the behaviour of the frequency 
starting off vertically upwards with pressure (see EIG.[^. 
The latter point implies that the pressure (force) may be 
inaccurately calculated under tensile stress (pressure re¬ 
mains at about -2GPa when the unit cell volume keeps 
increasing), where vdW plays an important part. We will 
resolve the former point and we avoid the latter point by 
focusing on the compressive part. 


C. Non-hydrostatic compression 
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FIG. 1. (colour online) [Modelling] Graphite under hydro¬ 
static pressure. The frequencies of the graphite Eiu and 
Ejg' are plotted against pressure. The pressures are plotted 
against the unit cell volumes at which they are calculated. Vb 
is the unit cell volume of unstrained graphite. The linear fit 
of the phonon frequencies at compressive pressure up to 10 
GPa is presented as dashed lines, black for Ei„ and blue for 


E 


( 2 ) 
29 ■ 


Modelling non-hydrostatic condition helps to investi¬ 
gate problems found under hydrostatic condition. The 
only study of graphite under non-hydrostatic condition 
reported so far is that of Abbasi-Perez et al.. They sug¬ 
gested that the contribution to the shift of the in-plane 
phonon frequency from the out-of-plane compression is 
so little that it can be neglectecpSl. The following results 
do not agree with that. 

We model uniaxial strain along the c-axis of graphite 
by varying the interlayer distance while fixing the in¬ 
plane geometry. The phonon frequencies and out-of- 
plane stress are calculated at each interlayer distance 
and fixed sp^ bond length. EIG. (a) shows the shifts 
of the in-plane phonon frequencies against out-of-plane 
stress along c-axis. In FIG.[^(b), the stress as a calcula¬ 
tion output, is plotted against the input — the interlayer 
distance in this case. The shift rates with stress up to 
about 10 GPa, by least square linear fits, are 0.8 and -0.2 
cm“^GPa“^ for Ei^ and E^g^ respectively. It is worth 
noticing that the elastic constant Cia, determining the 


Poisson’s ratio Vzx, is poorly defined due to the struc¬ 
tural anisostropy of graphite, but can be considered to 
be close to zercP^. Our calculated C 13 value is -10.5 GPa, 
and the corresponding is -0.024. Therefore, uniaxial 
compressive strain here induces in-plane tensile stress. 
The degeneracy of the two modes can be again seen in 

this case when the graphene sheets are pulled apart. The 

(2) 

problem is that the Ei^ and E^g' modes shift with oppo¬ 
site signs. 

Next we consider uniaxial stress on graphite along the 
c-axis, by varying the interlayer distance and optimiz¬ 
ing the in-plane geometry at each interlayer distance. 
The phonon frequencies, the out-of-plane stress and the 
sp^ bond length are calculated at each interlayer dis¬ 
tance. The effect of the negative Poisson’s ratio can now 
be clearly illustrated in FIG. as the in-plane bond is 
also compressed as we compress along the c-axis. The 
amount, however, is tiny. FIG. (c) presents the in¬ 
plane phonon frequency against out-of-plane stress and 
again the output stress is plotted against the input inter- 
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FIG. 2. (colour online) [Modelling] Graphite under uniaxial strain / stress along the c-axis. The frequencies of the graphite 
Ei„ and are plotted against out-of-plane stress in the case of uniaxial strain (a) and unaixial stress (c). The out-of-plane 
stresses are plotted against the interlayer distances at which they are calculated, in the case of uniaxial strain (b) and uniaxial 
stress (d). The linear fit of the phonon frequencies at compressive stress to about 10 GPa is presented as dashed lines, black 
for Ei„ and blue for Ejg^. 


layer distance in FIG. (d). The shift rates with stress 

up to 10 GPa in this case are 1.3 and 0.3 cm“^GPa“^ 
( 2 ) 

for Ei„ and , respectively. The shift rate with uni¬ 
axial stress for the Ei„ is about a quarter of the shift 
rate under hydrostatic stress; this is large enough to be 
significant. 

It is usual to present frequency against stress, be¬ 
cause stress is considered as input in experiments. How¬ 
ever, the atomic positions (strain) determine properties 
such as phonon frequency, and it is therefore helpful to 
plot the frequencies against strain. To be specific, for 
graphite, the shifts of the frequencies of the in-plane 
modes Ei^ and E^g are considered to be induced by in¬ 
plane strain. So we plot the phonon frequencies against 
in-plane strain under hydrostatic and non-hydrostatic 
conditions and then obtain the corresponding Griineisen 
parameters and 7^(2) for the Ei^i and Ef^ modes, 

respectively, according to Eq.|^ Values for the two modes 
should be the same from the two dimensional nature of 

Eq-lIHll 

In the case of uniaxial strain, the phonon frequencies 
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FIG. 3. (colour online) [Modelling] Graphite under uniaxial 
stress along c-axis. The sp^ bond length is plotted against 
the interlayer distance, at which it is calculated. 


for both modes shift at fixed in-plane geometry. The 
Griineisen parameters are hence 'Je, =00 and 7p,(2)=oo, 
according to Eq. ^ 
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In the case of uniaxial stress, in FIG. (a) we plot 
the in-plane phonon frequencies against the sp^ bond 
length, which is calculated by optimizing the in-plane 
geometry at each interlayer distance. The top axis of 
in-plane strain is converted from the sp^ bond length 
by £ = (a — ao)/ao x 100 %, where £ is the in-plane 
strain, a is the sp'^ bond length and og is the sp'^ bond 
length of unstrained graphite. This is the same data as 
in FIG. We apply a linear fit to the data points un¬ 
der compression up to about 10 GPa (the same as in 
FIG. [^, and obtain the shift rates dojE i^ /de=-152.00 
cm“^/% and du ( 2 )/de=-35.50 cm“^/%, which corre- 

^2g 

spond to 7 ei„= 4.86 and 7 ^( 2 ) =1.14. 

In the case of hydrostatic pressure, in FIG. H (b) we 
plot in-plane phonon frequency against the bond 
length, which is calculated by the geometry optimiza¬ 
tion at each unit cell volume. The top axis of in-plane 
strain is converted in the same way as before. The data 
is the same as in FIG.j^ We apply a linear ht to the data 
points under compression up to about 10 GPa (the same 
as in FIG. and get the shift rates du)/de=-00.20 
cm“^/% and doo ^ 2 )/de=-50.50 cm“^/%, corresponding 

to 7 b = 2.21 and 7 „( 2 )= 1 . 81 . 

^2g 

We present the Griineisen parameters obtained in the 
case of uniaxial stress and hydrostatic pressure in FIG. 
and compare them to that of uniaxial strain and those 
reported in TABLE The values for the Ei„ and E^^^ 
modes obtained under hydrostatic condition are in good 
agreement with the most values of the E 2 g mode of 
graphene presented in TABLE |l] Similar agreement has 
been reported in the literature; this is the quantitative 
reason why the in-plane phonon frequency shift s are con¬ 
sidered as induced by in-plane strain alon^^^^^^HIIl How¬ 
ever, the difference between the values of the Ei„ and 

f2') 

^ 2 g modes increases under uniaxial strain when the ra¬ 
tio of out-of-plane strain to in-plane is considerably larger 
than it is under hydrostatic pressure and neither of the 
values agrees with the E 2 g of graphene. Under uniaxial 
strain, where there is only out-of-plane strain, the differ¬ 
ence of the frequency shifts can be considered as infinity. 

The out-of-plane strain is responsible for the difference 

( 2 ) 

of the Griineisen parameters for the Ei^ and E^^ modes. 
Regarding the amount of the difference, to include the 
out-of-plane strain contribution to the in-plane phonon 
frequency is desirable in the case of hydrostatic pressure 
and definitely necessary in the cases of uniaxial strain 
and stress. 

To quantify this contribution, we refined Eq. ^as Eq.|^ 
The solutions to the secular equation of Eq. [6 led to the 
new parameter 7 ', as shown in Eq. for small shifts. 
Alternatively, we can rewrite Eq. [^with a full hydrostatic 
term {Sxx T ^yy 

“ ^(.^xx~^^yy~^^zz)^ ^yy) (T ^)^zz 

OJQ Z 

( 8 ) 

We believe this out-of-plane contribution is mostly re¬ 
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FIG. 4. (colour online) [Modelling] Graphite under uniaxial 
stress along c-axis (a) / hydrostatic pressure (b). The data 
in (a) is the same as in FIG. 2 (c) and the data in (b) is 
the same as in FIG. 1. Here the frequencies of the graphite 
Ei„ and are plotted against in-plane sp bond length, 
which is calculated at each interlayer distance (a) / unit cell 
volume (b). The top axis of in-plane strain is converted from 
the sp^ bond length by e = (a — ao)/fflo x 100%. The linear 
fit of the phonon frequencies at compressive stress to about 
10 GPa is presented as dashed lines, black for Eiu and blue 
for Ejg^. The Griineisen parameters for both modes, obtained 
from the linear fit, are presented in each case for comparison 
to the values in TABLE [I] 

fated to the compression of the tt electrons, which is be¬ 
yond the picture of the force constant model. 

Let us now return to Eq. In the case of uniax¬ 
ial strain, where the shifts of the frequencies are entirely 
from the out-of-plane strain, we plot the in-plane phonon 
frequencies against the interlayer distance — the calcu¬ 
lation input, in FIG. (a) and fit the data under com¬ 
pression up to about 10 GPa (the same as in FIG.[^(a)). 
The top axis of out-of-plane strain is converted from in¬ 
terlayer distance by £ = (033 — a 33 o/a 33 o) x 100 %, where 
£0 is the out-of-plane strain, 033 is the interlayer dis¬ 
tance and 0330 the value of unstrained graphite. The 
shift rates for the Ei„ and E^g modes are duJEi^lde=- 
0.915 cm“^/% and 9a; ( 2 )/3£=0.204 cm“^/%, corre- 
spending to 7 ^ =0.0585 and 7 ' ( 2 , =-0.0131, according 
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to Eq. They are small, but non-negligible as the out- 
of-plane strain is about 30 times larger than the in-plane 
strain under hydrostatic condition (from the anisotropy 
of graphite) and can be even larger under non-hydrostatic 
conditions. It is worth noticing that the in-plane phonon 
frequency cannot be considered as an indicator of the 
in-plane bond stiffness in this case as the Ei„ and E^^ 
modes, both representing the in-plane bond stiffness, 
shift with opposite signs under out-of-plane compressive 
strain. Now we have quantified the out-of-plane strain 
contribution by 7 ', which is responsible for the separat- 
ing of the Ei„ and E^g modes and then the in-plane 7 can 
be the same in various conditions for the two modes (and 
the E 2 g of graphene) as it should be from its definition. 

Finally, we calculate the refined value for the in-plane 
7 from Eq. 7 by 7 ^^ and 7 ' ( 2 ) obtained under uniaxial 

strain. For Doth modes, in the case of uniaxial stress 
and hydrostatic pressure, we calculate the contribution 
to the shifts of the frequencies from out-of-plane strain 
by y and attribute the rest to the in-plane strain. And 
from that we obtain the refined in-plane 7 . 

We plot the in-plane phonon frequencies against inter¬ 
layer distance, which is the calculation input under uni¬ 
axial stress in FIG. (b) and calculated by the geometry 
optimization at each unit cell volume under hydrostatic 
pressure in FIG. (c). The top axis of the out-of-plane 
strain is converted from the interlayer distance in the way 
as mentioned above. The data is the same as in FIG. [2] 
(c) and for uniaxial stress and hydrostatic pressure, 
respectively. We obtain the shift rates for the two modes 
at the same pressure range as before. Under uniaxial 
stress, we obtain 7 ^. =1.84 and 7 ( 2 ) =2.01 from the re- 

^ ^2g 

suits shown in FIG. (a) and[^ (b). Under hydrostatic 
pressure, the values are 7£:^,^=1.85 and 7 ^( 2 ) =1.90, from 

the results shown in FIG. (b) and[^ (c). 


D. Summary of the results and their applications 

The E 2 g and Fit, will be separately discussed so the 
subscript of the Griineisen parameter is removed in this 
paragraph. For E^g , the uniaxial strain modelling gave 
7 / = —0.0131. Applying this 7 / to uniaxial stress mod¬ 
elling we obtained 7 = 2.01 and to hydrostatic pressure 
modelling we obtained 7 = 1.90. The discrepancy is due 
to non-linear relationship shown in FIG.[^(a). We choose 
to use 7 = 1.90 because the data of hydrostatic pressure 
modelling was validated by comparing to experiments. 
Now Eq. can be used to give the value of 7 ° for each 
case in TABLE |T] with the approximation that graphene 
and graphite have the same elastic constants. The com¬ 
parison is presented in the table. No firm explanation 
can be made of the current results of graphene as the 
elastic constants (especially C 13 and C 33 ) of graphene 
cannot be accurately obtained and also due to the er- 
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FIG. 5. (colour online) [Modelling] Graphite under uniax¬ 
ial strain (a) / uniaxial stress (b) / hydrostatic pressure (c). 
The frequencies of the in-plane modes of the graphite Ei„ 
and E 2 g'^ are plotted against interlayer distance, which is the 
calculation input in the case of (a)&(b) and is calculated at 
each unit cell volume in the case of (c). The top axis of 
out-of-plane strain is converted from interlayer distance by 
£ = (<133 — 0330/®33o) ^ 100%. The linear fit of the phonon 
frequencies at compressive stress to about 10 GPa is presented 
as dashed lines, black for Ei„ and blue for Ejgh 


rors in these experiments, but it is clear that graphite 
has a smaller y than graphene because of 7 /. For the 
studies of graphite under uniaxial compression along c- 
axipISH, (e^a; + eyy)i the contribution of 7 / be¬ 

comes much more significant. Refinement can be done to 
include the mentioned non-linear effect but should await 
further work — either initial experimental evidence or, 
theoretically, the evolution of the 7r-electrons can be vi¬ 
sualised under compression to provide a clearer picture 



















of the phenomenon to be quantihed. For Ei^, the uni¬ 
axial strain modelling gave 7/=0.0585. Applying this 7 / 
to uniaxial stress modelling we obtained 7 = 1.84 and to 
hydrostatic pressure modelling we obtained 7 = 1.85 — 
excellent agreement achieved. The Griineisen parameter 
of Ei„ was considered the same as E^g and indeed they 
are close. But the 7 / of opposite sign for these two modes 
makes them distinguishable under hydrostatic pressure 
and further under uniaxial compression. 

Other analysis can also be done to study the phase 
transition induced by the change of interlayer distance, 
which has clear signs in the presented results, namely the 
significant drop of phonon frequency in FIG. i(b),|4](a) 
and (b) under large compression. 

IV. CONCLUSION 

We model uniaxial strain, stress along c-axis and hy¬ 
drostatic pressure on graphite and calculate the vibra¬ 
tional frequencies of the in-plane modes derived from the 


graphene ^ 2 g mode. The shifts of the frequencies come 
from both in-plane and out-of-plane compression. We 
quantify the contribution from out-of-plane strain by new 
parameters 7/_B2g=-0.0131 and 7/£;j^=0.0585, and there¬ 
fore refine the existing values of the Griineisen param¬ 
eter 7 as 7£;2g=l-90 and 7£;j^=1.85. This contribution 

( 2 ) 

is responsible for the separating shifts of the Ei„, E^g 
modes of graphite and the E 2 g of graphene under hy¬ 
drostatic pressure and therefore non-negligible, against 
previous conclusion. It can be significant under non¬ 
hydrostatic condition. A more reliable value of the in¬ 
plane Griineisen parameter is useful for strain calibra¬ 
tion in various applications and can be further refined by 
studying the 7r-electron behaviour. 
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